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Abstract 

We consider escape from chaotic maps through a subset of phase space, the hole. Escape rates are 
known to be locally constant functions of the hole position and size. In spite of this, for the doubling 
map we can extend the current best result for small holes, a linear dependence on hole size h, to in- 
clude a smooth h 2 lnh term and explicit fractal terms to h 2 and higher orders, confirmed by numerical 
simulations. For more general hole locations the asymptotic form depends on a dynamical Diophantine 
condition using periodic orbits ordered by stability. 

1 Introduction 

Here we consider open dynamical systems, in which motion is considered only until the trajectory reaches 
a specified subset of phase space, the "hole(s)". The initial conditions are distributed with respect to some 
measure, typically an invariant measure on the phase space of the corresponding closed system (the "escape" 
problem), or on such a measure restricted to hole(s) ("recurrence" or "transport" problems). There are also 
relevant questions where there are two systems connected by a hole (the "metastability" problem). There is 
a vast literature in mathematics and physics on open dynamical systems and applications, some of which is 
mentioned in gHEH2]. 

For the escape problem, we expect that for strongly chaotic systems (eg with exponential decay of 
correlations), the survival probability P(n) = fjt[M n ) decays exponentially with time n, so that an escape 
rate can be defined: 

7= lim lnP(n) (1) 

n— } oo fi 

Here fi is the measure of initial conditions (assumed to be invariant under the dynamics) and M n the subset 
of the phase space M that survives for at least n iterations before reaching the hole H C M. In general the 
limit might not exist, and we need upper and lower escape rates, possibly infinite. Here we restrict to maps, 
considering only a discrete time n; similar questions apply to flows. 

It is of interest to consider how the escape rate depends on the hole, in particular its size h = fJ,(H). On 
one hand, there have been attempts to quantify the small hole asymptotics of the escape rate, that is, its 
behaviour for a sequence of holes contracting to a point x £ M. Ref. [T7] numerically observed the effects 
of periodic orbits on the escape rate from a chaotic map with a small hole. Ref. [5] expresses escape rates 
in terms of correlation functions, finding good numerical agreement in the case of diamond billiard with 
continuous time escape rate computed to order h 2 ; the heuristic arguments proposed there suggest that for 
billiards with holes typically considered (holes that allow particles to escape that reach a small region of the 
boundary, but from any direction) the first order result is impervious to a short periodic orbit located in the 
hole, but that corrections would be needed for the second order formulas. Ref. [Hj shows rigorously that 

1 = h(l-A- 1 )+o(h) (2) 

for a wide class of hyperbolic maps, where A x is the stability eigenvalue of the periodic orbit starting from 
x € M, infinite if the orbit of x is aperiodic. This result (amongst others) was also shown independently in 
a study of the open doubling map [6j. 

On the other hand it is known that the largest invariant set that never reaches the hole, a subset of M^, 
is a locally constant function of the hole. This is also mentioned in [TJ] but was noted much earlier [3T]. It 
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then implies the escape rate, and also properties such as the relevant Lyapunov exponents and entropy, are 
locally constant. This can be seen (but was not remarked on) for the diamond billiard in Fig. 4 of Ref. [S]. 
While putting an apparent limit on small hole asymptotic results, this effect is also related to the very helpful 
observation that open dynamics restricted to the repeller is normally a subshift of finite type; almost always 
in the case of the doubling map considered here. For very recent results on the existence, continuity and 
locally constant properties of the escape rate for more general hyperbolic systems, see Ref. [7J. 

We note that the locally constant property does not exclude the possibility that the escape rate is smoother 
as h — > than the single derivative required by Eq. pi). This paper is an attempt to push the small hole 
asymptotics as far as possible in the simplest possible dynamical system, the doubling map, with the hope of 
elucidating similar effects in more general contexts. Sec. [2] discusses previous work relevant to the doubling 
map. Sec. [3] considers a hole starting at zero and uncovers an explicit, but fractal, expansion for the escape 
rate to arbitrary order in h, including a smooth h 2 In h term as the first contribution beyond linear. Sec. [4] 
considers arbitrary hole locations, finding numerical support for a general h? In h term for rational locations, 
but rigorous arguments for terms as slow as h/lnh for some Liouville locations, ie very well approximablc 
by rationals. In more general maps, the expansion appears to depend on an approximation theory using 
periodic orbits ordered by stability. Finally, Sec. [5] explores some connections with the Gauss map. 

2 Existing results 

The doubling map x —¥ 2x (mod 1) has an invariant measure /i simply given by Lebesgue, so the hole size h 
is simply the length of the relevant interval(s). There is a correspondence between the binary representation 
of a point x and symbolic dynamics for the partition {[0, 1/2), [1/2, 1)}, modulo minor details to do with the 
dyadic rationals (ie fractions with denominator a power of two). 

The case of a hole between and h has received some attention in the literature, also including more 
general expanding circle maps [131 121] . In brief, given a hole size h, find 

h- = sup{2"x > x (mod 1), Vn G N} (3) 

x<h 

This is closely related to the concept of a "minimal prefix" [16] . The quantity h- is almost always a dyadic 
rational i/2 n for some integer i, and in this case h lies in the interval of constant escape rate i/2 n < h < 
i I (2™ — 1) , often abbreviated here as just "interval." Note that each periodic orbit of the map appears exactly 
once as the upper limit of one of these intervals. 

Within an interval, the dynamics is a subshift of finite type. If we write the (n term) binary expansion 
of 1 - h: 

n 

l-/^=$> fe 2- fc (4) 
fe=l 

then the leading eigenvalue of the transition matrix is /3/2 where f3 is the solution of 

n 

l = $>fc/3- fe (5) 
fe=i 

with )3«2, Thus we find the escape rate 

7 = -lnf (6) 

Both the local dimension of the set of irrational /i_ and the Hausdorff dimension of the relevant repelling 
set are equal to In/?/ In 2 = 1 — 7/ In 2. For irrational h- the escape rate is computed as in the rational 
case (but taking n — > 00) since it is continuous and monotonic. Thus the case of H — [0, h] is of interest for 
/3-expansions. 

3 Small hole asymptotics 

We now consider the question as to what asymptotic results are possible in this case. The intervals of 
constant escape rate give an important constraint on the possible form of such asympototics: 7 is constant, 
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but h has varied by an amount proportional to h 2 (roughly equal to h 2 when i = 1). Thus jh~ 2 has jumps 
of unit size occurring over a short distance relative to h itself: the 0(h 2 ) term in the expansion, if it exists, 
must have jump discontinuities. 

For any h contained in an interval, we have j(h) = j(h-) so it is sufficient to consider just values h- on 
the left of an interval. From the monotonicity and continuity we may furthermore calculate 7 as a limit of 
rational or irrational h- as convenient. For the following argument we use irrational h = h— and the binary 
expansion of h rather than the 1 — h above: 

00 

h = 2- 1 K a = 2- , Y J A k 2~ k (7) 

k=l 

Here, J = (In K — In K)j In 2 is a non-negative integer so that K < 1. There are at least two logical choices 
for J: We can set 1/2 < Kq < 1 so that Kq (and each K m below) is bounded above and below; this is most 
helpful in determining the magnitude of the relevant terms. The other choice is to set J — which simplifies 
the computations. For now we need to keep J arbitrary. The other new quantities in the above expression 
are A k = 1 — o-j+k, with all the non-positive A k equal to zero. Thus the expansion for /? becomes 

00 00 

l = Y f a k p- k = (J3-l)- 1 -Y f A k r J - k (8) 

k=l k=l 

We know from [14] that to leading order 7 h/2, since the fixed point at x = has A = 2. Thus we 
write 

00 

= 2- 5^ (9) 
1=1 

where ei = 1, and we allow the other e; to be functions o(/i _1 ) as h — > 0. Substituting for /3, expanding both 
powers of (3 in a binomial expansion, substituting h/K for 2~ J and equating coefficients with equal powers 
of h, we find that 

<> = 1 ~, T7~~ (1<» 

1 ( K 2 2K\ llXi 
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and so on, where 

00 

K m = Y, A * km2 ~ k ( 12 ) 

k=l 

consistent with the definition of Kq above. The K m for m > are fractal functions bounded on [1/2, 1) as 
shown in Fig. 1; they are discontinuous at dyadic rationals. 

We also see that at second order the expansion involves Jh 2 which gives rise to a smooth h 2 In h term. 
Thus we can extend Ref. [T3] result to 

h h 2 \\nh\ h 2 ( 3 \nK KA 
^2 + ^lnV + T(-2 + lnf + ^ 

for h an irrational value in the fractal set not contained in any interval. 

If h lies in an interval [h-,h+] — [i/2 n ,i/(2 n — 1)] we identify it with the discontinuity occurring at 
Ka = i/2( n ~ J \ The jump at this point is — = 2^ 1- ^™ - ' 7 ^ independent of i, which leads to a 
cancellation to second order with the discrepancy in the first term, since /i_ — h+ = i/2 2n to this order; all 
other terms are smooth. Thus, the jumps in 7 are exactly those required by the locally constant intervals. 

In what sense does this expression converge in the h — > limit? Considering a fixed domain of Kq, 
intervals containing 2~ J Kq decrease with increasing J. Thus we expect convergence for almost all fixed 
Kq G (1/2, 1) as J — > 00, where h — 2~ J Ko, with an error of smaller order than the last retained term. On 
the other hand, there is no guarantee that for fixed J the terms with higher powers of h continue to decrease. 

Recall that we did not specify J exactly; we can see that the transformation J — > J + 1 has the effects 
Jf — > Kq/2 and K\ — > (Ki + Kq)/2, leaving the above expression for 7 unaffected. But this also means 
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Figure 1: The K m functions defined in Eq. (12 1 



that the coefficient of h 2 is log-periodic, ie a periodic function of In h. We can thus write the expansion of 
7 as a series involving h, h 2 \nh and /i 2 + 27rl1 ?/ ln 2 for q e Z, however the coefficients (Fourier expansion of 
the K functions in the variable ln h) do not appear to have a simple analytic form. They are also slow to 
converge due to the discontinuities. We can continue analogously with ft, 3 and higher, indicating that the 
Mellin transform of 7 has a semi-infinite periodic forest of poles reminiscent of those for Julia sets [3] . 

Using the arbitrariness in J, we now set J = to yield the simplest form of these expressions: we can 
substitute h for K and evaluate the other K functions at h, giving the main result: 



7 
7i 

72 
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71 + 72 + 73 + 0(h 4 In 3 h) 

h 

2 

2K x h - 3/i 2 
8 

6Kf h + 3K 2 h 2 - 21K 1 h 2 + Uh 3 



(14) 
(15) 

(16) 

(17) 



where we note that K m (h) ~ h\ In™ h\/\a m 2 as h -> 0. 

Fig. [2] shows a comparison between the escape rate calculated directly and with the above formulas. 
The direct numerical computation uses a Markov matrix of side length N — 2 16 (N — 2 19 for h < 0.01) 
that does not need to be stored due to its regular structure. The leading eigenvalue corresponding to 7 is 
obtained by repeated multiplication of the matrix on a vector with initially identical nonzero entries, using 
long double (19 digit) precision. The direct escape rate is plotted, then subtracting each of the ji in turn, 
leading to smaller and smaller remainders. We can see that the locally constant intervals (smooth curves in 
the direct computation) become smaller for smaller h, but remain as spikes at dyadic rationals. Elsewhere, 
the remainder decreases with the order of the approximation, showing good convergence despite the fractal 
character of the functions involved. 

Various connections may be made between the K functions and sclf-affinc functions and sets appearing 
in the literature. The graph of K\ is easily seen to be a discontinuous self-affine set, see Ref. [19j for some 
relevant properties. It shares the relation 



K x {x/2) = (K 1 (x)+x)/2 



(18) 



with the Takagi function, used to study diffusion in extended doubling maps |15) , but the similarity appears 
to end there, as the Takagi function is continuous. We can construct a generating function for all the K 
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Figure 2: Numerical test of the small hole expansion for the escape rate, Eq. ( 14) 



functions, (cf Eq. [5]) 



OC 

k=l 



(19) 



so that, at least formally, we have 

K n (x) = (-0dp) n G (x)\^ 2 (20) 

The limits here are quite subtle, since Gp(x) effectively computes the Bernoulli convolution corresponding 
to /3, and it is known that the distribution of the latter is singular when (3 is Pisot, including a sequence 
accumulating at f3 = 2 20J; it might be interesting to investigate these connections further. 



4 Dependence on hole position 

For holes at general locations, H — [Hi,H2] with H2 = Hi + h we no longer have an explicit formula for 
the characteristic polynomial (except for Markov holes; see [11]) so it is not yet possible to carry out the 
same expansion as above. However there are a few observations that can be made, which also apply to more 
general piecewise expanding maps. 



Heuristic argument for h 2 In h There are reasons to believe the h 2 In h term is more widely applicable. 
In a general hyperbolic map, the linear dependence of the escape rate on the periodic orbit stability arises 
because that is the extent of the "shadow" [I] cast by the hole on itself under iteration around the periodic 
orbit. Beyond this, we need to study the shadow of the hole on itself along orbits that do not remain close 
to the short periodic orbit (if there is one) passing through the hole. The shortest of these are typically of 
length |lnft|, with stabilities of order h^ 1 . This will make a piece of the hole of size roughly h 2 inaccessible 
to the escaping trajectories. At the next length there will be more periodic orbits, balanced by greater 
instabilities, also yielding of order h 2 . This continues until double the original length, say 2|ln/i| beyond 
which the periodic orbits are effectively just shadowing combinations of shorter periodic orbits. From this 
argument we might expect these to give a total contribution to the escape rate of order h 2 In h. However it 
will depend on the details of the periodic orbits. We now introduce a more systematic approach and show 
that reality is more subtle. 



Relation to dynamical Diophantine properties Let us fix Hi and take h — > 0. In general we expect 
that the escape rate is given (roughly) by 

7^1-A- 1 ) (21) 
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where A is the stability eigenvalue of the shortest periodic orbit in the interval. The approximation arises 
from omitting the small h limit needed for Eq. If Hi is itself aperiodic, we have y/h — > 1 as h — > 0, 
however a short periodic orbit x lying just above h\ will lead to smaller 7 when h > x — Hi; if there are 
infinitely many such orbits (where "short" is relative to an appropriate function of h), this may affect the 
limiting escape behaviour beyond linear order. 

So, assume that H\ is not itself periodic and consider a sequence x n of periodic points approximating H\ 
from above. We obtain a sequence 

7 „«/ l „(l-A- 1 ) (22) 

of approximate escape rates with hole sizes h n — x n — Hi. The correction to the linear term, /i^A" 1 is of an 
order of magnitude relative to h n that depends on the stability of the periodic orbit. 

(Pre-)periodic hole locations For a rational Hi = p/q, that is, a periodic or preperiodic points for the 
doubling map, the distance to a rational approximation (ie periodic point) with denominator 2" — 1 cannot 
be smaller than l/(q(2 n ~ 1)), thus the size of a hole covering Hi and its approximation is bounded by a 
known constant multiple of A -1 and so the correction term (making the same convergence assumption) is 
order h 2 , hence not dominating the heuristic h 2 hih term above. A symbolic version of the same argument 
could also extend to more general expanding maps. It may be possible to find explicit formulas for the 
characteristic polynomial in simple cases (piecewise linear map, short (pre)-periodic orbit), leading to the 
coefficient of h 2 In h along the lines of the previous sections, however we will confine ourselves here to a 
numerical investigation of the doubling map. 

We consider hole locations Hi which are all the rational values with denominators in {1, 2, . . . , 10, 15, 21, 31, 63} 
thus including all periodic points with periods p < 6 and also 13 preimages of orbits with periods p < 4. The 
matrix size used was N = 9999360 = 2 10 x 3 2 x 5 x 7 x 31 and the power method continued until a tolerance 
of 3 x 10 -14 was reached, using 19 digit precision. For each location, the escape rate was computed for holes 
of size h — 41/iV and h = 2 P 41/N, corresponding to each other under the p times iterated map. The number 
41 was chosen as having no common factor with the denominators, and roughly the smallest size consistent 
with the above tolerance. Removing the known linear term from the escape rates, these two points are then 
used to estimate h 2 In h and h 2 terms corresponding to the sequence of hole sizes found by multiplying by 
powers of 2 P . As with the hole at zero discussed previously, the coefficient of the h 2 Inh term is found to be 
approximately independent of h, but the coefficient of h 2 is a fractal function of h. 

The results, depicted in Fig. [3j are consistent with the formula 

7 = (l-A- 1 )fe+ (1 ~ A ^" h 2 \\nh\+0(h 2 ) (23) 
In 2 

where A is taken to be infinite for preperiodic points, and the In 2 is chosen for consistency with the previous 
case Hi = 0. The single outlier near the top right of the figure corresponds to the point Hi = 62/63, the 
period 6 orbit most affected by the fixed point at 1; here convergence of the escape rate to the small hole limit 
is slowest, so the extrapolation is not as accurate. In more general cases, note that where a periodic orbit 
has negative stability eigenvalue (as in the tent map for odd periods) Eq. ([2]), and hence its generalisation, 
need to be modified if the periodic point lies on the boundary of the hole. 

Typical hole locations Typical behaviour is given by Borel-Cantelli results, which as discussed in Ref. pQ 
apply to the doubling map, and more generally. The relevant theorem states that an infinite number of the 
iterates 2 n x (mod 1) of the point x under the doubling map are contained in a shrinking sequence of 
intervals I n for almost all x if and only if the sum of the measures (ie lengths) of the intervals diverges. 
For example, the intervals [Hi — (n In 7i) -1 , Hi] contain infinitely many y n — 2 n Hi (mod 1) for almost 
every Hi if and only if k < 1. Each such point leads to a periodic point x n = Hi + (hi — y n )/(2 n — 1) 
lying in [Hi, Hi + ((2™ — l)nln n.) -1 ]. Thus we expect an infinite sequence of holes h n of size as low as 
((2™ — l)nlnn) _1 with escape rate 

7« ~ M 1 - 2 ~") ~ M 1 - h n \lnh n \\n\\nh n \) (24) 

This argument is not rigorous as it depends crucially on the convergence rate of Eq. ([2]) but suggests that 
for typical Hi a bound on the correction term might need to be slightly greater than \hr \nh\. 
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Figure 3: Numerical investigation of the dependence of the ft 2 |ln/i| term for rational hole locations; see 



Eq. (23) 



The well approximable case In contrast to the above heuristic and numerical arguments, the following 
negative result is rigorous: 

Theorem 1. Given any continuous strictly increasing function g defined on [0,1/3] with g(0) = and 
g(l/3) > 1/4, there is an irrational (ie aperiodic) Hi such that the doubling map with hole [Hi, Hi + h] has 
escape rate j(h) so that "f(h)/h — 1 approaches zero but is not 0(g(h)) as h — > 0. 



Proof: Consider 



x k 



* 1 
l-^ 2 nk - 1 



(25) 



fc=i 



for a strictly increasing sequence with ni = 2. Clearly all the xk lie in the unit interval. Then, we define 



Hi = lim xk 



(26) 



and construct the sequence so that Hi satisfies the required property. 

First, assume that divides Uk+i for all k. Then, 2" fc — 1 divides 2™ fe + 1 — 1 and so xk is a rational with 
denominator 2™ A ' — 1 (or a factor) and hence a periodic orbit of length uk (or a factor). Furthermore, since 
the sequence is strictly increasing, n^+i > 2nk which implies that the limit Hi is irrational as required. 

There are two further lower bounds on the growth rate of n^. One comes from the limit in Eq. ([2]): There 
exists a function H(e) (implicitly depending on xk) such that for all h < H(e) we have 



|7(J0/ft-(l-2-*)|<e 



(27) 



considering the hole [xk — h, x K\ shrinking to the periodic orbit xk- We choose e = 2~ 2 ™ fc to ensure that 
7(/i)//i is sufficiently close to that determined from the periodic orbit. This requires nk+i > — log 2 H(2~ 2rik )+ 
1, where the final term ensures that the hole can extend downward beyond xk+i to h\. 

Finally, we need to ensure that the correction 2~" fc is greater than g(h) for the hole [hi,xx]- This is 
accomplished by requiring nk+i > — log 2 G~ 1 (2~" fc ) + 1 where G _1 is the inverse of G, an arbitrarily chosen 
monotonic function with the same conditions as g but with G(h) ^ 0(g(h)) as h — > 0. Thus any sequence 
satisfying divides rifc+i and 



>ik- 



-i > max[— log 2 H(2 



~2n k 



l,-log 2 G-\2- n *) + l] 



(28) 



satisfies the requirements of the theorem. 



□ 
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Examples: For g(h) = h a with < a < 1 we need at least rik+i > nk/a. For g(h) — — 1/ln/i (so that 
j(h) = h + 0(hj In h)) we need at least nk+i > C2 nk , which has extremely rapid growth leading to Hi as a 
Liouville number. 

Remarks: The main idea of the above proof is to find aperiodic points well approximated by periodic 
orbits of the map, in particular such that the separation is much smaller than the inverse stability of the 
periodic orbit, A -1 . Approximation by periodic orbits has also appeared in different contexts, for example 
in obtaining the multifractal properties of Lyapunov exponents [18j . Ordering periodic orbits by stability is 
a more generally useful technique for systems with widely varying expansion rates |10j . 



5 Epilogue: The Gauss map 

Finally, there is an intriguing connection with h 2 \nh in the literature, namely for the Gauss map x —> {!/%} 
with invariant measure dfi — (In 2(1 + x))~ 1 dx. The natural symbolic dynamics x —> |1/ X J gives the 
partial quotients of the continued fraction expansion x = [0; a\a2 ■ ■ ■}■ The periodic orbits are quadratic 
irrationals. Thus, a hole H = [0, l/(n + l)] corresponds to escape when > n. The hole size is h = J H d[i = 
ln(l + (n + l) _1 )/ln2. Hcnsley [H] analyses this problem using the transfer operator 

giving its leading eigenvalue A(s, n) in his Eq. (7.9) as 

/ 9 \ /n 2 6- 12\ 1 6hxn 

A ( 2 ~n' n J =1+ (^2hi2~ J n ~ ^2 + ^ 

Setting this equal to unity, he finds that the value of 9, and hence the Hausdorff dimension of the invariant 
set, to contain a term proportional to lnn/rt 2 , ie also h 2 \ah. The escape rate, given by — lnA(2,n), has 
however no such term! 

For more general hole positions, the previous discussion suggests that the asymptotics depends on the 
approximation of the hole position using periodic orbits ordered by stability. The Gauss map has infinitely 
many branches, however there are only finitely many periodic orbits with a given stability bound. It would 
thus be interesting to analyse the escape rate for general hole positions using approximation by quadratic 
irrationals. 
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